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3 ' Abstract 

O^, We study the effect of small decoherence in continuous-time quantum walks 

on long-range interacting cycles (LRICs), which are constructed by connect- 
CN ■ ing all the two nodes of distance m on the cycle graph. In our investigation, 

each node is continuously monitored by an individual point contact (PC), 
1/-^ I which induces the decoherence process. We obtain the analytical probability 

On ' distribution and the mixing time upper bound. Our results show that, for 

small rates of decoherence, the mixing time upper bound is independent of 
I distance parameter m and is proportional to inverse of decoherence rate. 
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1 Introduction 



Random walks on graphs have broad applications in various fields of math- 
ematics, computer science and natural sciences, such as mathematical mod- 
eling of physical systems and simulated annealing [1]. The quantum me- 
chanical analog of the random walks on complex networks has been studied 
with respect to the localization and delocalization transition in the presence 
of site disorder [^lElll]. Quantum walks (QWs) have been largely divided 
into two standard variants, the discrete-time QWs (DTQWs) O [6] and the 
continuous-time QWs (CTQWs) [7j. The DTQWs have been investigated 
on trees [8], on random environments [9], for single and entangled parti- 
cles [lOj and also in [TTl [12]. In the recent years, the CTQWs have been 
studied on n-cube [13], star graph OdS], small- world network [16], quo- 
tient graph [17] , line [El [I9l [20] , dendrimer [21] , distance regular graph [22] , 
circulant Bunkbeds [23], odd graph [21] and on decision tree [25l [26], 

In all of these cited works have been supposed that we have a closed quan- 
tum system without any interaction with its environment. Firstly, Kendon 
and Tregenna considered the effect of decoherence in quantum walks in " De- 
coherence can be useful in quantum walks" [27]. By numerical observation, 
they found that a small amount of decoherence can be useful to decrease 
the mixing time of discrete quantum walks on cycles. Then deoherence in 
quantum walks has been studied on line [28], on circulant [29] and on hyper- 
cube [30j . In " Mixing and decoherence in continuous-time quantum walks 
on Cycles" |31] . it has been provided an analytical counterpart to Kendon 
and Tregenna's result for the continuous-time quantum walk on cycles. Its 
results show that, for small rates of decoherense, the mixing time decreases 
linearly with decoherence while for large rates of decoherence, the mixing 
time increases linearly toward the classical limit. Moreover, for the middle 
region of decoherence rates, the numerical data confirms the existence of a 
unique optimal rate in which the mixing time is minimal. 

In this paper, we consider the effect of decoherence in continuous-time 
quantum walks on long-range interacting cycles (LRICs) as the extensions 
of the cycle graphs. LRICs are constructed by connecting all the two nodes 
of distance m on the cycle graph (nearest-neighboring lattice). A detailed 
description of the network structure will be given in the next section. Nu- 
merical analysis of CTQWs on LRICs has been provided in "Coherent exci- 
ton transport and trapping on long-range interacting cycles" [3]. We take a 
small amount of decoherence into account and by analytical technique show 
it can decrease the mixing time in continuous-time quantum walks on 
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Figure 1: Long-range interaction cycles G(8;3) and G(10;4) 



LRICs. For this end, we use of Gurvitz model ^32j. In this model, 
each of vertex is monitored by a corresponding point contact induced the 
decoherence process. We calculate the analytical probability distribution 
then obtain the mixing time upper bound for small rates of decoherence. We 
show that it is independent of the distance parameter m and is proportional 
to inverse of decoherence rate. Our paper is structured as follows: After 
a detailed description of the network structure LRICs in Sec. 2, we study 
continuous-time quantum walks over the underlying structures in Sec. 3. 
In Sec. 4 we consider the effect of decoherence in CTQWs on LRICs and 
in Sec. 5 we focus on the decoherent CTQWs when the decoherent rate is 
small. In Sec. 6, we define the mixing time and obtain its upper bound. 
Finally, in Sec. 7 the conclusions are presented. 

2 Structure of LRICs 

To construct long-range interacting cycles (LRICs), we can use of the fol- 
lowing rule: First, we construct a cycle graph of N nodes where each node 
connected to its two nearest neighbor nodes. Second, all the two nodes of 
distance m on the cycle graph are connected by additional bonds. LRICs are 
denoted by G{N, m) that N is the network size and m is long-range interac- 
tion parameter. LRIC is a one-dimensional lattice with periodic boundary 
conditions and all nodes of the networks have four bonds [3j. The structures 
of G(8,3) and G(10,4) are illustrated in Fig. 1. 

3 CTQWs on LRICs 

CTQWs on LRICs is obtained by replacing the Hamiltonian of the system 
by the classical transfer matrix, i.e., H = —T [25j. The transfer matrix T 
relates to the Laplace matrix by T = — 7A, where for simplicity we assume 
the transmission rates 7 of all bonds to be equal and set 7 = 1 in the 
following. In Laplace matrix A, nondiagonal elements Aij equal to 1 if 
nodes i and j be connected and otherwise. The diagonal elements Ajj 
follow as Ajj = —ki that ki is the degree of vertex i. The basis vectors |j) 
associated with the nodes j of the graph span the whole accessible Hilbert 
space. 
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Then, the Hamiltonian matrix H of G{N, m) {m > 2) takes the following 
form, 



' -4, ifi=j; 
1, i{i = j± 1; 



1, ii i = j i: m; 
0, Otherwise. 

The Hamiltonian acting on the state \ j) can be written as 



(1) 



H\j) = -4|i) + \j - 1) + I J + 1) + \j - m) + \j + m) , (2) 

which is the discrete version of the Hamiltonian for a free particle moving 
on a lattice. It is well known in solid state physics that the solutions of the 
Schrodinger equation for a particle moving freely in regular potential are 
Bloch function [331 [M]. Thus, the time independent Schrodinger equation 
is given by 



H|$e) = Eel^e), (3) 

where the eigenstates |<I>e) are Bloch states and can be written as a linear 
combination of states |j) localized at nodes j, 



The projection on the state \j) is $6)(j) = {j\^e) = which is 

nothing but the Bloch relation ^e{j + 1) = e~*^^>0(j) [331IM] . 
Now the energy is obtained from Eqs. (3) and (4) as 

= -4 + 2cos(6l) + 2cos(m6'), (5) 

for j = 0,1,..., N — 1. The classical and quantum transition probabilities 
between two nodes can be written as 



PkAt) = J2e-'^'{km{^e\j), (6) 
e 

7:,,{t) = \akAt)\' = \Y^e-^''^o{k\^e)mj)\'. (7) 

e 
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Figure 2: Fig. 2(a) shows point contact detector j monitoring the electron 
in dot j and Fig. 2(b) shows point contact j in the presence of the electron 
in dot i + 1. 



4 The Decoherent CTQWs on LRICs 

Gurvitz's model 

To analyze the decoherent continuous-time quantum walks on LRICs, we 
use analytical model developed by Gurvitz [32l [35| . 

In this model, every node is regarded as a quantum dot. Thus LRIC is 
represented by a set of the identical tunnel-coupled quantum dots (QDs). 
The walks are done by an electron initially placed in one of dots. A ballistic 
point-contact (PC) is placed near every dot that is taken as a noninva- 
sive detector. We assume that all point contacts are identical. Also, they 
are placed far enough from quantum dots so that the tunneling between 
them is negligible. Moreover, for simplicity, we consider electrons as spin- 
less fermions [31j . Each PC continuously monitors the attached quantum 
dot. This measurement process induces decoherence to electron walks as 
is shown in Fig. 2. Firstly, we study the simple quantum walks with the 
Hamiltonian [36] 

(8) 

^ -^A,,(t)(|f)(i| + U)(^))+^6,(t)|j)(j|, 

«i i 

where |j) = c^|0) denotes the state that the electron is placed at dot j. 
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Aij is the hopping amphtude between dots i,j and ej{t) is the on-site dot 
energy. We assume the constant hopping amphtude between hnked dots and 
no on-site terms. For simphcity, we renormahze the time, so that it becomes 
dimensionless |37] . Hence, the Hamiltonian has the form 

1 

-^0 = 4 5Z + 4^i+l + + (9) 

j=0 

The tunnehng Hamiltonian Hpc j describing electron transport in the 
point contact j can be written as [371 132] 

^PC,j = X] J + ^ ^rjal jCirj + X flirj (ajj^/ J + H.C.). (10) 

/ r Ir 

Here, and aj! j(ar,j) are the creation (annihilation) operators in the 

left and right reservoirs of detector j, respectively, and Qirj is the hopping 
amplitude between the states I and r of detector j. 

The presence of an electron in the left dot results in an effective increase of 
the point-contact barrier (fi^^ — > Qi^^. + (Jfi^ ,,), we can represent the interac- 
tion term as [37l [32] 

Hint,j = X dftirjc'jCjialjarj + H.C.). (11) 
Thus, the entire system Hamiltonian can be described by 

N~l 

H = Hq+ X Hpcj + Hintj- (12) 

j=0 

We suppose for simplicity that the hopping amplitude of j-th point- 
contact is weakly dependent on the states / and r, so that it can be replaced 
by its average value, {i^irj — f^) and Silirj — 5^1. The occupation of the 
quantum dot can be measured through the variation of the detector current 
AI = I2 — Ii where Ii = elnQ'^ pijprjVj is the detector current when the 
electron occupies the first dot, Fig. 2(a) and I2 = e27r{Q'^ + 6Cl)'^ pijprjVj 
is the current flowing through the detector in the presence of the electron 
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in the second dot, Fig. 2(b). The density of states in the left and right 
reservoirs of detector j are pij and prj and the voltage bias is the variation 
of the chemical potentials in the left and right reservoirs in detector i.e. 
Vj = fiij — iJ,r,j [35j. Now, we are ready to write Schrodinger equation for 
the entire system with Hamiltonian H. The effect of the detector on the 
quantum dot can be obtained by tracing out the detector states. Gurvitz has 
shown, for a double-dot and detector together (Fig. 2), Schrodinger equation 
results in the evolution of reduced density matrix traced over all states of 
detector which coincides with the Bloch-type rate equations [32l[35]. These 
equations are as 

= ^^o(Pi+i,i - Pj,j+i), (13) 

Pjd+i = ^^jPj,j+i + ^^o{Pj,j - Pj+i,j+i) - ^Pj,j+i, 

where €j = Ej — Ej^i and 0,q is the coupling between the left and right 
dots. Also, pjj+i(t) = p*^i j{t) are the off-diagonal reduced density matrix 
elements and the diagonal terms of this density matrix pjj{t), pj^ij^i(t) 
are the probabilities of finding the electron in jth. dot and in j + 1th dot, 

respectively. Moreover, T = {\J^ — \f^Y^ is decoherence rate due to 
continuous observation with a non-invasive detector [32l [35] . Applying this 
model to our system results in 



^Pj,k{i) — f [— /5j-l,A: — Pj+l,k — Pj-m,k — Pj+m,k + Pj,k~l 



(14) 



+Pj,k+1 + Pj,k-m + Pj,k+m] - r(l - 5j^k)Pj,k 

Our subsequent analysis will focus on the variable Sj^k defined as 

Sj,k = i^'^PjM (15) 
Substituting Eq. (15) into Eq. (14), we have 

do 1 r Q _i_ Q A— m+1 Q Am+1 a c 

'dt ~ iL '^j-l,k ~r i^j+l,k ' >-'j-m,k * >~'j+m,k ^j,k-l 

(16) 

+Sj,k+i + i'^^^Sj^k-m + - r(i - 5j,k)Sj,k 
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5 Small Decoherence 



We consider the coherent continuous-time quantum walks when TN <C 1. 
Eq. (16) can rewrite as the perturbed hnear operator equation [31] 

±S{t) = {L + U)S{t), (17) 
where the hnear operators L and U are 

(18) 



= - Sa,p)Sa,,S^,u, (19) 

that L is a A^^ x N"^ matrix and is the entry of L indexed by the row 

index (/i, i^) and the column index (a,/3). Also, U has the same behavior. 
By the above Substituting, we obtain 

p^,,= (20) 

where < a, f3, fijU < N — 1. The initial conditions are 

Pa,f3{0) = 5'a,/3(0) = Safl^/Sfi. (21) 

To obtain zero-order solution of Eq. (20), we require to an expansion on 
the eigenvectors of L^^^'^l or 

N-l 

where < A;, / < iV — 1. After some algebra with Eq. (22), one can obtain 

■K{k + l) TT{k-l) TTm{k + l) /^"^(^ - M /oo\ 

Afc,/ = i[sm{^j^) cos( ' ^ ] + « [sm( ^) cos( '-)] (23) 

and 

. 7r(A^+/) Ti{k-l) .^+1 . 7rm(fc + . 7rm(fc - /) 
Afc,/ = 2[sm( — — )cos( — — — + [sm( — )sm( — )],(24) 
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when m is an odd and even number, respectively. Eigenvectors of L are 
given by 

Now, we consider eigenvalues of the unperturbed linear operator L (Eqs. 
(23), (24)) carefully. We investigate the important degeneracies of the eigen- 
values Xk^i of L that lead to non-zero off-diagonal contribution of U. Firstly, 
due to symmetry of Eq. (23), we have A^^; = Xi^k, while their eigenvectors 
(Eq. (25)) are clearly different (Note that in Eq. (24) there is not such 
symmetry). The second subset of degenerate eigenvalues appear when we 
replace k + l = 0{modN) in Eqs. (23) , (24). one can see the corresponding 
eigenvectors are not the same. The third subset of the degenerate eigenval- 
ues reveal when we set k = I. 

First-order correction terms are given by the diagonal elements of Eq. (19) 
calculated on eigenvectors of Eq. (25). For the first subset of degenerate 
eigenvalues, they are equal to — F^^^^y^. By introducing this eigenvalue per- 
turbation to each pair of A^j,^;) with k ^ I and k + l ^ N, the degeneracy 
of the first subset is removed. The degeneracy of the second subset is ab- 
sent in our case since the their eigenvectors are anyway excluded from the 
final solution by initial condition [37J. In our discursion, there is not the 
degeneracy of the third subset since U is diagonal over the corresponding 
eigenvectors [3T]. For these eigenvalues, the correction terms are —T ^^^^^ . 
Thus, as mentioned in (SUET], Eq. (20) can be expressed in terms of eigen- 
vectors of Eq.(25) 

c, 5a,0 1 ST t-i X X \ 27ri(A;a + 1(3) 
Sa,/3[0) = ^ + ]^ (1 ~ ^fc+'.o ~ ^k+i,N) exp[ — ] (26) 

k,l=0 

The full solution is of the form 

k,l=0 

The probability distribution of the continuous-time quantum walks is given 
by the diagonal elements of the reduced density matrix Pj{t) = Sjj{t), that 
is 



N-l 



k,l=0 



'Jk+1,0 - 
iV2 



-ik+LN 



exp[ 



j2iri{k + l)j 



]e*(-^(fc,o+^{fc.o) (28) 
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Now, we want to rewritten the above equation for odd and even values 
of m, respectively, 
odd m 

The solution is of Eq. (20) is 



PA^) = F + 2^ ]y2 ^ * + (1 - 4,/)e ^ iv *]e iv 

/c,/=0 



X exp[it sin(^) cos(^) + i^t sini^^^) sin(^Ilf^)]. 

for odd m. 
even m 

The full solution of Eq. (20) is 



fc,'=o (30) 
X exp[itsin(^^) cos(^^) + i^+H sin(I!!lMO ) sin(II!l(JzO )]. 
for even m. 



6 Mixing time 

To define the mixing time of continuous-time quantum walks, we use the 
principal motivation to studying random walks. In computer science, the 
probabilistic algorithm provide the best solution for many problems. Thus 
the precise solution is obtained by a well-chosen sampling distribution. Gen- 
erating the such distribution is a matter of mapping the uniform distribution 
into the desired one. Hence, it is important to get a truly uniform distri- 
bution [38] . The mixing tiniG T^nix 

(e) is defined as the number of steps 
needed before the distribution is guaranteed to be e-close to the uniform 
distribution, or 

N-l ^ 

T^Ue) = min{T : ^ |P^.(t) - _| < e}, (31) 

3=0 

that Pj{t) is the probability distribution of quantum walk on node j of 
graph G and jf is the uniform distribution over graph G. 
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Based on the above analysis, we will be interested in decreasing the mixing 
time. 

Firstly, we want to calculate the upper bound on the e-uniform mixing time 
T^Ue) of Eq. (29) (for odd m). 
We define 



^^■(*) = iv ^ exp[^^sm(— ) +z™tsm(^^)]e— . (32) 

fc=0 



Hence, we have 

Af-l 



Af-l N~l 

2TTi(k + l)j ^ ^ -I — ^ 2-Kik{2j) 



k,l=0 k=0 

Hence, Eq. (29) can be rewritten as 



m - ^1 = e-^^*|M|(-) - - + ^^[M,,{t) ]|.(33) 

Since |Mj (t)| < 1, we obtain 



(34) 

< e-r^*|l + ^(l-|)|. 
Based on the definition of time mixing. 



Yl - < e-^^*(iV + - 1) (35) 
i=o 

_rt p iv-2 . 

Because of e ^ < 1 , the above equation reduces to Ne n < e. 

As a result we obtain the mixing time upper bound of 

IN 2 
T™.(e)<-ln(-)[l + ^^]. (36) 

This relation is in agreement with the result is mentioned in |31j for cycle. 
Note that in this case, the mixing time lower bound can not derive of Eq. 
(33) easily, since there is a relation between M?(|) and M2j{t). 
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Now, we calculate the upper bound on Tmixi^) of Eq. (30) (for even m). 
We define 

Mj{t) = -^Yl exp[i-(sin(— ) + sin(— ))]e^^ 

k,i=o (37) 

X exp[-i™+if (cos(^) - cos(^))]. 

Thus, we have 

\m)-^\= e-rV*|M,(t)-^| 

(38) 

< e-rV*[l + ^] 

and with summation over j and using the mixing time definition, we get 
Ne'^—^ < e. 

This gives the mixing time upper bound of 

1 1 iV 

TmiAe)<f[l + j^]H-)- (39) 

The lower bound of mixing time can derived by the first equality of Eq. 

(38), as 

\Pjit)-jr\= e-rV*|M,(t)-i| 

(40) 

= 

where in the last inequality, we set Mj{t) = In other words, in this time 
the quantum walk completely reaches to uniform distribution 
Thus, the lower bound of mixing time follows as 



Tmixie) > 0. (41) 

Eqs. (36) and (39) show that the Tmix upper bound is independent of 
the distance parameter m. Moreover, since we approximated the coefficient 
of exponentially function in Eqs. (34), (38), the mixing time upper bound is 
exactly proportional to ^ that accord with Fedichkin, Solenov and Tamon 
's result for the continuous-time quantum walks on cycles [31]. Also, these 
relations prove that the mixing time upper bound for even m is smaller than 
mixing time upper bound for odd m. 
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7 Conclusion 



We have studied continuous-time quantum walks on long-range interacting 
cycles (LRICs) under small-decoherence TN <C 1. We obtained the proba- 
bility distribution analytically and found that the mixing time upper bound 
for odd values of m {Tmixi^) < y l^(T)[Ar^]) larger than the mixing time 
upper bound for even m {Tmix{^) < r I^It)!!^^!])- These relations show 
that the Tmix upper bound is inversely proportional to decoherence rate T 
and is independent of the distance parameter m. Also, we proved that for 
even m the lower bound time mixing is zero. In other words, we have shown 
that Tmixif) decrease with F at least as fast as l/F. 
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